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by 

Ashok  K.  Kalra,  James  A.  Showalter  and  Lisa  A.  Pflug 
Naval  Research  Laboratory 
Stennis  Space  Center,  MS  39529-5004 


ABSTRACT 

We  have  processed  the  archived  basebanded  con^lex  time  series  hydrophone  data,  over  the 
sediment  pond,  that  was  acquired  with  high  resolution  waveform  Upsweep,  HFM,  210  - 
280  Hz).  This  data  was  first  unbaseband^  to  its  original  band,  i.e.,  from  210  to  280  Hz 
and  then  matched  filtered  with  the  source  replica.  For  the  purpose  of  geological 
interpretation,  it  was  then  basebanded  to  a  lower  frequency  which  improved  the  time 
resolution  of  the  data.  The  coherent  low  velocity  noise,  that  was  present  throughout  the 
data  set,  was  later  removed  with  proper  velocity  filt^.  These  processing  stqis  helped  in 
revealing  the  vertical  incident  coherent  events  and  we  give  our  interpretatic»i  of  these 
events. 

INTRODUCTION 

We  have  processed  the  high  resolution  ^psweep,  HFM  210  -280  Hz)  hydrophone  data 
(not  beam  formed  data)  of  Track  5,  that  is  located  over  the  sediment  pcmd,  in  order  to 
facilitate  the  geological  interpretation  of  the  sediment  pond.  In  all.  Track  5  consists  of  30 
pings,  from  181  to  210,  with  a  ping  interval  of  1582  m  but  only  selected  pings  (188  to 
192, 195, 196, 201, 204  to  208)  were  available  for  processing.  Although  128  channels 
(including  two  desensitized  phones)  of  hydrophone  data  were  aa)uiral  for  each  ping,  for 
some  pings  only  64  channels  were  available  for  frirther  processing  and  analysis. 

PROCESSING  OF  HYDROPHONE  DATA 

The  hydrophone  data  available  on  tape  consisted  of  complex  time  series  which  had  been 
complex  demodulated  to  save  storage  space.  Consequently,  the  first  task  was  to 
unbaseband  the  complex  time  series  data  to  its  origimd  frequencies  (210  -  280  Hz)  and 
convert  it  to  the  real  hydrophone  data. 

The  complex  time  series  ^ta  was  read  using  the  NRL/SSC  Off  Line  Processor  (OLP)  and 
converted  to  SEGY  (Society  of  Exploration  Geophysicists  'Y'  format)  using  UNIX 
equipment  and  modified  software  developed  by  Scri]ms  Institution  of  Oceanography.  This 
data  was  still  in  complex  demodulated  form,  llie  SEGY  data  was  written  to  exabyte  tapes 
for  storage.  Additional  software  was  written  by  Richard  R.  Slater  (NRL^SC)  to  convert 
the  SEGY  format  to  ASCII  for  portability.  This  data  which  was  still  in  the  complex 
demodulated  format  could  now  be  read,  unbasebanded  and  converted  to  real  data  by 
software  written  in  PV_WAVE  environment  (Pflug,  Caruthers  and  Slater,  1993). 

It  was  decided  that  the  first  twenty  seconds  of  hydrophone  data  was  more  than  sufficient 
for  later  (geological)  interpretation.  The  original  ctnnplex  basebanded  data  was  sampled  at 
128  Hz  and  so  2560  sample  points  corresponding  to  20  seconds  were  further  processed. 

To  return  the  complex  demodulated  data  to  its  original  real,  high  frequency  btmd,  the  data 
was  shifted  up  in  frequency  (frequency  modulated)  to  the  original  210  -  280  Hz  region  and 


paired  with  the  corresponding  negative  frequency  complex  conjugate  mirrcn’  image.  This 
required  that  the  sampling  rate  be  increased  to  a  minimum  of  560  Hz.  The  modulated  real 
data  was  interpolated  to  a  sampling  rate  of  640  (or  12  8(X)  teal  sample  points  for  twenty 
seconds),  a  rate  for  which  the  data  is  somewhat  oversampled. 

The  real  data  was  matched  filtered  using  the  original  210  -  280  Hz  (Upsweep,  HFM) 
source  replica.  The  time  zero  was  defined  to  be  the  (xigin  of  source  transmission. 
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Figure  1 


Figure  2 


At  this  stage,  a  display  of  the  high  frequency  wideband  (210  -  280  Hz)  data  indicated  that 
time  domain  interpretation  of  such  a  high  fi^uency  data  was  difficult  (see  Figure  1).  In 
order  to  improve  ^e  interpretational  aspect  of  the  data,  the  real  data  was  base^ded  to 
lower  frequency  of  0  Hz.  The  additional  step  of  basebanding  the  data  involved  shifting  the 
210  -  280  Hz  positive  and  negative  frequency  regions  down  to  0  Hz  and  applying  an 
inverse  Fourier  transform  to  obtain  time-domain  data.  This  resulted  in  a  simpler  time 
domain  wavelet  and  increased  the  time  resolution  of  the  data  set  (Figure  2). 

DISCUSSION  ON  PROCESSING 

Basebanding  the  data  to  lower  fr^uency  results  in  a  relatively  higher  rado  of  main-tobe  to 
side-lobe  energy  in  time  domain  in  comparison  to  the  one  basebanded  to  higher  frequency. 
This  is  illustrated  in  Figure  3.  Figure  3  shows  two  time-domain-wavelets;  one  basebanded 
to  0  Hz  (5  -  75  Hz)  and  another  at  205  Hz  (210  •  280  Hz)  with  an  identical  bandwidth  of 
70  Hz.  The  wavelet  basebanded  to  higher  frequency  (Figure  3b)  exhibits  a  number  of  side 
lobes  that  are  quite  high  in  amplitude  as  compared  to  the  main  lobe.  The  wavelet 
basebanded  to  lower  frequency  (Figure  3a),  on  the  other  hand,  indicates  only  one  side 
lobe  on  either  side  of  the  main  lobe  and  the  side  lobe  energy  is  much  lower  than  the  main 
lobe.  This  contributes  to  the  ease  in  the  interpretational  aspect  of  the  wavelet  basebanded  to 
lower  frequency  such  as  the  one  in  Figure  3a. 
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Figure  3 
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Figure  4  Figure  5 


This  was  further  verified  by  a  wedge  model  which  consists  of  two  geological  horizons, 
forming  a  wedge,  indicating  a  varying  degree  of  separation  (thickn^)  between  the 
horizons  from  one  channel  to  the  next  Figures  4  and  5  are  the  time  responses  to 
unbasebanded  (210  -  280  Hz)  wavelet  and  waveltt  basebanded  to  0  Hz  (S  >  75  Hz) 
reflectively.  CHjviously,  the  response  to  wavelet  basebanded  to  lower  frequency  O^igure 
5)  is  better  able  to  detect  the  existence  of  two  closely  spaced  horizons  and  thus  exhibits 
better  time  resolution. 


Next,  data  basebanded  to  0  Hz  was  velocity  filteied  to  enhance  the  vertically  arriving 
events.  Figure  6  shows  the  data  of  Figure  2  after  velocity  filtering,  and  it  clearly  has 
enhanced  the  vertically  arriving  events  and  eliminated  the  interfering  coherent  evrats 
arriving  along  the  hydrophone  array. 


INTERPRETATION  OF  DATA 


Processed  data  (such  as  of  Figure  6)  indicates  a  lot  of  wedge-type  horizons  interfering  with 
(Hie  another.  However,  we  do  not  expect  such  smictures  to  be  real  based  on  geolo^cal 
inferences  of  the  depositional  environment  in  this  area.  This  is  an  area  where  pelagic 


HVtMOmONE  DATA  FOR  TRACK  S 
■ASeBANOED:S-TS»« 
VELOCITT  filtered 
PINa  No.  Ill 


T—wIMt  — — »  ■  If  H  a  M  u  n  u  ••  la  in 


Figure  6 


sediments  accumulate  at  rates  measured  in  centimeters  per  million  years.  These  are 
fine-grained  sediments  (mostly  clay  with  some  minor  carbonate)  derived  from  the  settling 
of  particles  from  the  water  column.  Consequently,  we  do  not  expect  any  significant 
acoustic  impedance  contrasts  within  the  sediments  to  generate  the  reflection  events.  We 
give  below  our  interpretation  of  the  data. 

This  data  was  acquired  with  a  towed  vertical  array  source,  consisting  of  10  elements  with  a 
separation  of  3.66  m  and  the  hOTizontal  towed  receiver  array  with  a  near  offset  of  912  m 
and  far  offset  of  1224.5  m  (ARSRP,  1991).  With  the  water  depth  of  4130  m  from  the 
towed  arrays  and  for  the  flat  sea  floor  case,  the  incident  angles  at  the  receiver  array  from  6  ** 
to  8**  approximately.  The  radiation  pattern  from  the  source  within  these  angles  vary 
considerably  especially  considering  the  presence  of  source  and  receiver  ghosts.  In  practice, 
there  was  also  coupling  between  the  each  source  elements  which  further  complicated  the 
signal  radiated  by  the  source  in  the  vertical  direction.  It  was  not  possible  to  deghost  the 
data  due  to  the  complicated  source  radiation  pattern.  Consequently,  we  believe  that  the  data 
of  Figure  6  represents  in  interference  pattern  of  the  complicated  source  signal  (including  the 
ghosts)  reflected  back  from  the  sea  floor. 

CONCLUSIONS 

It  is  illustrated  that  the  time-resolution  for  the  purposes  of  geological  interpretation  of  the 
hydrophone  data  is  improved  when  the  data  is  processed  with  wavelet  for  which  the 
fi^uency  band  is  shifted  to  start  at  0  Hz  (basebanded  to  0  Hz)  in  comparison  to  the  one  for 
which  the  passband  starts  at  higher  h-equency.  This  should  be  true  for  beam  forming 
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